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1 .  Introduction  and  St 


Several  types  of  applications  entail  point  processes  with  "dead" 
tines  after  each  event.  One  such  fanily  of  stochastic  processes  is 
fl(PRP),  the  fanily  of  Pareto  Renewal  Processes.  The  i.i.d.  inter¬ 
arrival  tines  (Xj)  satisfy,  F(x)  ■  P(X  <  x)  »  1  -  (£)  ,  x>A>0 
and  s  >  0. 

An  additional  interesting  property  of  the  interarrival-tine  distri¬ 
butions  is  that  they  are  all  "thick- tailed"  relative  to  the  corresponding 
distributions  for  Gaussian  processes  and  Poisson  processes.  Further, 
a  variety  of  tail  thicknesses,  one  for  each  s-value,  is  possible. 


These  two  properties  lead  to  sone  interesting  inference  problem,  of 
which  one  is  here  concerned  only  with  signal  detection. 

[The  Pareto  distribution  itself  was,  of  course,  introduced  by 
Vilfredo  Pareto  (1848  -  1923).  (See  Reference  [22]).  This  distribution 
has  been  used  and  studied  by  mnwrous  other  authors  including  Pigon  (1932) 
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Mandelbrot  (1960,  1963),  Fisk  (1961)  end  Johnson  and  Kotz  (1970).] 

Both  one~saaple  and  two-saaple  signal  detection  problems  with 
historical  data  will  be  considered  here.  The  organization  of  the  paper 
is  as  follows.  Section  2  contains  the  basic  properties  of  the  inter' 
arrival-tine  distributions  and  the  MLE's  (next was  likelihood  estiaates) 
imder  various  circuastances.  Section  3  contains  the  distribution 
theory  necessary  for  insight  and  inference.  Moot  of  the  proofs  for 
the  results  of  this  section  are  straightforward.  In  Section  4,  one 
introduces  the  fundanental  statistical  concepts  to  be  used,  naaely, 

(i)  the  BDT  (Basic  Data  Trans  tarnation) ;  (ii)  the  M-S-S  (aininal 
sufficient  statistic);  (ill)  the  M-S-N  (naxiaal  statistical  noise); 
(iv)  POP  (parasMtric  distribution- free)  statistics;  and  (v)  MFDF 
(nonparaswtric  distribution- free)  statistics. 

Section  5  discusses  the  uses  of  the  Kolaogorov-Saimov  (1933) 
statistic;  and  its  extensions  by  Ulliefors  (1967,  1969),  Srinivasan 
(1970)  and  Y.  Choi  (1980)  in  signal  detection.  Sections  6  and  7  treat 
the  one-sanple  and  two-saaple  detection  p rob lens ,  respectively.  Tables 
suanarizing  the  results,  are  presented  at  the  end  of  appropriate  sections. 

Appendix  A  contains  proofs  of  sons  of  the  assertions  nade  in  the 
paper.  Appendix  B  contains  mate ri cal  exaaples  illustrating  each  of 
the  detection  procedures  derived  in  Sections  6  and  7. 
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I 

2.  Ejnwtfflf  Progrtlw  and  Estina tors  of  Pareto  gwiwl  Proces s 

Most  of  tho  results  of  this  *  set  ion  are  adapts t  loot  of  results  in 
Johnson  and  Kotx  (1970). 

•ach  A  >  0,  s  >  0  by 

I 

F(*)  - 

^  This  is  a  special  fora  of  the  Pearson  typo  VI  distribution.  Tho  Pareto 

density  function  is 

f(x)  •  sA*  X'*'1  for  a  >  A  >  0. 

If  X  %  Pa(A,s),  then  tho  rth  nonont  of  X  exists  if  and  only  if 
r  <  s  and  is  given  by 

I  ^  • 

The  variance  of  X  is  sA*(s  -  l)"*(s  -  2)“l  for  s  >  2.  For  further 
I  details  on  nonents,  see  Malik  (1966). 

Malik  (1966,  1967)  has  also  obtained  results  on  the  characteristic 
functions  of  the  order  statistics  frou  a  Pareto  distribution  as  well  as 
I  recurrence  relations  between  the  aonents  and  covariance  of  the  order 

statistics.  Levy  (192S)  discovered  a  class  of  stable  laws  (Stable 
Pares i an)  which  follow  the  asynptotic  fore  of  the  Pareto  law. 

I 

» 


The  Pareto  distribution  Pa(A,s)  is  defined  for 


{ 


1  -  (7)  *  >  A 

0  x  <  A 


If  X  •v  Pa(A,s),  Y  •  In  X,  then  Y  'v  T-exp(ln  A,s)  «diere 
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T-exp(B.s)  if  truncated-exponential  with  distribution  taction 
H(*)  •  1  *  •“***  '  **  for  x  >  I.  Also  if  l  *  »*',  thon 
Z  \  11(0,  A*9).  Inference  for  thn  trwcated- exponential  has  boon  studiod 
by  Pori  end  U.  Oioi  (1978),  log  (1981)  and  signal  detection  p rob loo  for 
the  uni  fore  process  has  boon  investigated  by  Y.  Choi -Bell -Ahead -Part 
(1982).  Sow  of  those  results  will  bo  used  and  colored  to  the  ones  herein. 

Here  ow  is  primarily  interested  in  MLB's  of  the  perawtert  A 
and  s.  Other  types  of  estimates  are  given  in  Johnson  and  Kott  (1970). 

The  likelihood  taction  for  the  Interarrival  tlws  (Xj,  Xj . X||) 

froe  a  PUP  is  L  •  A  — 9*  ■  .  The  proofs  of  the  following  theorens 

J-l  t*  1 

are  straight  forward. 

Theorea  2.1.  (One-parawter,  •  *  *0  known) 

(i)  The  MU  of  A  is  A  •  X(l),  which  is  distributed  Pa(A,  nsQ); 

(ii)  E(A)  -  nSgA(ns0  -  l)"1  for  os0  >  1  and  V(A)  .  nSgA^ntg  -  2)(n*0  -  l)2)*1 
for  n*0  >  2; 

(ill)  The  MVUE  (ainiflw  variance  unbiased  estiwtor)  of  A  is 
A*  •  (ns0  -  l)(ns0)“,X(l)  fOr  ns0  >  i,  and  V(A*)  «  (ns0(ns0  -  2)]'lA2 
for  nsQ  >  2, 

(iv)  Both  A  and  A#  are  consistent  estiaators  of  A,  i.e.,  A/A  *  1, 

A*/A  ♦  1  as  a  ♦  •  but  only  A*  is  unbiased. 

(v)  the  M-S-S  (ainiaal  sufficient  statistic)  for  A  is  X(l)  and 
the  faadly  (Pa(A,as0))  indexed  by  a  >  0  is  complete. 


Theorea  2.2.  (One -par water,  A  •  Ag  known) 
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(i)  Ttii  MU  of  i  is  s  • 


f  In  X.  •  n  In  A- 

J.|  >  " 


<“»  '  4..  *«i>  *  £  **•  *<7>  '  A  i 

»  .  »  »  n»* 


.V 


.....  «(.*.  •  <*).  '«»>>-  »(•>  * 


(iv)  Tho  NVUE  of  s  is  |ivtn  by  s*  •  (““Jo  ud  V(s*)  •  ~  ■ 1  j 
for  n  >  2; 

(v)  loth  s  and  s*  art  cons  is  tent  tstiastors  of  s  but  only  s* 
is  tssbiasod;  sad 

(vi)  tbs  M-S-S  for  s  is  s  {or  *•). 

Thoor—  2.3.  (TWo  parnatttrs ,  A.  s  both  unknown) 

(1)  Tho  MU  of  (A.s)  is  (X,s)  whom 


X  •  X(l)  and  s  - 


f  In  X.  - 

...  ' 


nln  X(l) 


(ii)  X 


s  art 


indopoodont »  X  ^  Pa(A,ns) ,  —  <v  x? 


2(n*l)  * 


(ill)  E(X)  •  -~~r  for  ns  >  I,  V(X) 


nsA 


(ns -2)  (ns-1) 


y-  for  ns  >  2; 


2.2 


(W)  E(s)  •  ^2-y  for  n  >  2,  V(s)  •  — 


0»*2)  (n-5) 


for  n  >  3; 


(r)  Tht  MVUE  of  A  is  A*  •  ^  X(l)  for  ns  >  1  and 

ns 

A2 

ns (ns -2) 


V(A#)  • 


for  ns  >  2; 


(vi)  til*  MVUE  of  s  i«  »* 


and  VC**) 
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.  h  -  vi. 

n 

(vil)  th*  M-S-S  for  (A,*)  i*  (i,s). 


* 


Th*  results  of  this  section  sr*  iwsriud  In  Table  2.1  below.  The  last 
colun  of  this  table  yields  th*  H-S-N  (naxinal  statistical  noise)  for 
various  detection  problems .  A  fornal  definition  of  th*  H-S-N  is  |iv*n 
in  Section  4.  As  far  as  Table  2.1  is  concerned,  on*  should  view  the 
M-S-N,  N(£),  as  coapleuentary  to  and  statistically  independent  of  the 
M-S-S,  S(^).  Several  versions  of  th*  M-S-N  are  given  In  Section  3. 


3.  Distribution  Theory  for  Pareto  general  Processes 

In  this  section  on*  develops  son*  results  which  yield  the  M-S-S 
(to  be  defined  in  Section  4)  for  various  one-sanpl*  and  two-saapl* 
problean. 

Let  Xr  Xj . \  be  i.i.d.  Pa(A,s)  and  let  3(1)  «  X(2)  <  ...  < 

X(n)  denote  the  order  statistics  of  th*  X*s. 

Definition  3.1.  On*  denotes  by  G-0-S(k)  th*  distribution  of  th* 
order  statistics  Induced  by  a  randoo  ssnple  of  site  k  drawn  fro*  a 
distribution  G(-).  The  following  leona  is  fondanental. 

1 flu  3.1.  (Mvi  (19S3) ,  Pyhe  (196S)).  Ut  C,.  C2 . Ky  be 


t».(,  },  lo.),  I*. I  »»*  *r)M  ia  V*t  lea  J,  kfitidai  I.;. 


i.i.d.  Exp(X).  Then  (i)  k£  »  \  r(k,X),  (ii)  £*  ■ 

j-1  5  * 

k-1 

C  C  ^ 

{  -i-  ,  -2 . -ill - }  %  U-0-S(k  -  1)  where  U(-)  is  the  U(0,I) 

kl  kl  k£ 

distribution ,  (iii)  k£  end  ere  indeoendent. 

Renarfc.  Note  that  X  has  been  eliminated  in  £*. 


3.2.  Let  Hj,  n2 . be  Exp(X)  -  0  -  S(k).  Define 


°! . "i  *  "f  •  ‘V  "i  *  "i  ♦  »  •  ‘'"i . "i  •  ni  *  •••  * 

nj_,  ♦  Ck  ♦  1  -  J)nj . ng  ■  kn  and  n##  • 


n!  ni  nf  . 

I  -  -  |  lhmtk 

t  ^  - .  1  *  ll,fn 


(i)  ng  %  r(k,X)  la  the  N-S-S  for  X; 

(ii)  n##  %  U-0-S(k  -  1); 

% 

(iii)  ng  and  n*#  are  independent. 


Theorem  3.3.  Let  Xj , 

X  JL 

(i)  (in  -i.  In  ,  ... 

(ii)  (In  X(l),  In  X(2) ,  .. 

(iii)  (in  In 


Xg,  X  be  i.i.d.  Pa(A,s),  then 
Xn 

,  In  --  )  i.i.d.  Exp(s) ; 

.,  In  X(n))  *v  T  *  E*p(ln  A,s)  -  0  -  S(n); 


Exp(s) 


0  *  S(n) ; 


•  •  •  •  9 


(**)  *x7nT  *  xfi-i iT  *  ••••  xHT*  *  •*•*»  ?*(*)  *•  th* 

to  x  <  o 

x*  o  <  x  <  1, 

1  X  >  1 

I  •  I  •  I  f  t 

(*)  U  1  •  Iffsrfyl  »  •••*  CjfyYyl  )  *  u(o,A”*)-o»s(n), 

(vl)  {im]  *  lIRTIT1  •  ••••  l^?TT1  ]  %  u-°-sW- 

Paflnltion  3.3.  Ut  Xj,  Xj . XQ  bo  i.i.d.  Pa(A,s).  (!*• 

following  randoa  variables  T.  9  .  f  ,  1  <  r  <  n  -  1  win  bo  uood 

T  r  T  mm  mm 

ft roqoontly  cho  root  of  tbo  paper  and  no  nontlonod  In 

Tablo  2.1.] 

(i)  T  •  f 

r  J«1 

(H)  Pr  -  )  ♦  («  -  r)  In  ,  1  <  r  <  n  •  1; 

(ill)  if  .  (n  ♦  1  -  j)  1  «  r  <  n  -  1. 

Th*  following  tbooron  is  an  easy  consequence  of  Lmms  5.1  and  5.2. 

Ti  r7  T.  » 

Theoron  5.4.  (i)  T  •  fyi*  .  .  yHlL  )  %  U-O-Sfn  -  I), 

<v  a  n  n 


1  <  r  <  »; 


and  is  independent  of  T  %  r(n,s); 


H  U-0-S(n  -  2)  and  is 


1 

I 

I 

I 

■* 

ti 

I 

I 

■ 

4 

. 

I 

i 

i 

\ 


T- 


.  JO  - 


independent  of  P  j  ^  I*(n  -  l.a); 

E.  f2  C».2 

(ill)  f  •  {  ,  -*= -  .....  p2-i  )  ^  U-0-S(n  *  2)  and  is 

'  i»-i  n-2  n-l 


independent  of  £  j  ^  l*(n  -  l.a). 


The  following  theores  contains  various  vara  ions  of  tha  N-S-M  (saa 
Da  fi  nit  ion  4.1)  dim  A  la  unknown.  Thf*  will  ba  us  ad  in  tha  Kolnoaorov- 
Sadmov  statistics  for  tha  one-saaple  casa  in  Saction  6. 


Thooren  S.S.  (i)  (In  («$&)*  . x  &*(•)• 


0-S(n 


-  I).  (II)  (ilfi  ,  lift . ito  1  v  •  1): 


(111)  (iff4}]’.  IiJ^TT)* . iffjfl*'  v  »-«*«•  ■  "• 

«»)  ufHb<<"n\  .  I4felil*>  «•.« i 


(»)  (a  la  iffy}).  (a  -  l)la  |fff}| . la  1.1.4.  e»(«). 


Da  fl  nit  Ion  3.3.  for  1  <  )  <  a  -  1,  ona  dafinas 

(1)  »,  •  lal^^l.  (ID  »,  •  ")  *  'tTS^T1’  • 

('  *,  *  llprl'"  «*>  2,  •  (»  *  I  '  ))  l«  t4^)- 

They  will  ba  oaod  throughout  tha  rest  of  the  oapar. 


11  - 


The  final  leaaa  for  this  section  deals  with  a  two -sample  situation. 


3.6.  Let  XJ#  X2 . Xm,  Yj ,  Y2,  . . . ,  Yn  be  independent  with 


X.  -v.  PafAj.Sj)  and  Y.  ‘v-  Pa(A2,s2).  Let  nL  ■  l 


X(j*l) 


n2  •  j  In^ffl  *  then  (*)  XO)»  Y(l),  Hj,  n2  are  independent; 


(ii)  X(l)  -w  Pa(Aj ,  asp,  Y(l)  ^  PafAj.nSj);  (iii)  Hj  <v  T(»  -l.Sj), 


n2  -v  T(n  -  l,s2);  (iv) 


In  and  2ns 2  In  -ffi-  are  i.i.d. 

A- 


1  A. 


X*.fcq«l/2);  (,)  2»,«,  '  »<•  Vl  *  2,2"2  ' 


X2||_4  where  N  ■  a  «  n;  and 


CN  -  2)  lea.  In  -  n*2  In 


2.^,  ♦  2s2n2 


^  F(2,  2N-4) . 


Beaarfc .  LesMa  3.6  (ri)  is  useful  when  a  ■  n  because  if  A}  •  A2 
end  $j  •  t^l*n  ***•  expression  on  the  left  does  not  involve  A  nor  s. 

(See  Section  7.) 

The  results  of  this  section  are  suanarized  in  Table  3.1  below. 

4.  Basic  Statistical  Concepts 


(A)  The  BOT  (Basic  Data  Trans foraatlon)  and  MSN 


Leal  Statistical 


Noise) . 


TABLE  3.1.  Distribution  Theory  for  Pareto  Renewal  Processes 


Table  3.1.  (Continued) 


-  is  - 


Let  the  generic  data  point  be  denoted  by  Z  *  (X. ,  . . . ,  X  )  in  the 

%  *  n 

one-sample  cases  and  by  Z  •  (Xj,  ....  X^,  Y^,  ....  Y^)  »  (Zj,  . ..,  Z^) 

in  the  two-sample  cases;  and  let  S(Z)  denote  the  M-S-S. 

Definition  4.1.  Let  W(Z)  be  a  (vector-valued)  statistic  independ- 
-  ^ 

ent  of  S(Z)  and  such  that  5(Z)  *  [S(Z),  W(Z)]  is  1-1  a.e.  Then, 

t\j  'V* 

(i)  5(0  is  called  the  BDT;  and  (ii)  N(Z)  is  called  the  M-S-N. 

It  is  known  that  S(Z)  contains  all  relevant  information  about  the 

*\j 

parameter  (vector);  and  it  will  be  seen  that  N(Z)  contains  all  relevant 

*\j 

information  about  the  structure  of  the  process.  From  5(£)  »  [S(Z) , 

N(Z) ] ,  one  should  almost. always  be  able  to  reconstruct  the  original 

data,  Z. 

«v. 

Example  4.1.  (i)  In  Lemma  3.1,  S(5)  *  kT;  W(£)  «  £*,  and 

“  'it  'Vi  'Xt 

6(E)  »  [kC,  5*]  is  the  BDT.  From  6(5),  one  can  reconstruct  5. 

'V  %  %  'w 

(ii)  In  Lemma  3.2,  S^J)  «  n*  and  W^n)  «  n**.  Hence  6(n)  *  (n£,  jQ**). 
From  the  examples  above,  it  should  be  clear  that  there  are  several 

possible  versions  of  the  H-S-S  and  M-S-N. 

Remark .  The  importance  of  these  concepts  is  that  as  a  rule  of  thumb 
in  a  goodness -of- fit  problem,  the  decision  rule  should  involve 
only  the  M-S-S  while  in  a  class-fit  problem,  the  decision  rule  should 
only  involve  M-S-N.  This  will  be  seen  in  the  sequel. 


(B)  Distribution-free-ness . 


There  are  two  types  of  distribution  free  statistics  that  arise 
in  detection  procedures. 


Definition  4.2.  (i)  A  statistic  T(Z)  is  called  nonparanetric 

distribution  free  (NPDF)  with  respect  to  a  family  ft*  of  stochastic 
process  laws  if  there  exists  a  single  distribution  Q(»)  such  that 

P(T(Z)  <  t | L>  -  Q(t)  for  all  Left*, 

'Xj 

(ii)  A  family  of  statistics,  (Tt(Z;  L)}  indexed  by  the  members 
L  e  ft*  is  called  parametric  distribution  free  (PDF)  with  respect  to 
ft*  if  there  exists  a  distribution  function  such  that, 

P{T.(Z;  L)  <  t)  -  Q.(t)  for  all  L  e  ft*. 

1  1 

It  is  clear  that  each  NPDF  statistic  is  PDF. 


1.2.  Let  ft*  *  ft(PRP) , 


first  N  interarrival  times  with  Z. 


and  Z  ■  (Zj ,  . . . ,  Zjj)  be  the 
■  Xj  for  1  £  j  £  m,  and 


Z  *  Y  for  1  <  r  <  n,  where  N  ■  m  ♦  n.  Let  s.  and  s_  be 

m  ♦  r  r  —  —  *  i  2 

the  MLE's  given  by  Theorem  2.2;  Tj  ■  — and  T2  ■  .  Then,  Tj 


is  NPDF  wrt  ft*  with  Q  *  F(2m,  2n) ;  and  T2  is  PDF  wrt  ft* 

2 

with  Qj  ■  X^.  Furthermore,  it  can  be  shown  that  Tj  is  a  function 
of  the  Z  only  through  the  M-S-N,  N(Z) ;  and  T_  is  a  function  of 
X  *  (X,,  ...,  X )  only  through  the  M-S-S,  S(X). 

1  m  r\j 

The  M-S-N  for  the  respective  cases  are  given  in  Table  2.1.  The 


relations  between  the  M-S-S,  M-S-N  and  DF  statistics  is  best  given 
by  the  following  theorem. 


Theorem  4.1.  Let  fl*  be  a  family  of  cdfs  admitting  a  M-S-S,  S(Z), 

■  *\» 

for  data  Z  «  (Z, ,  . . . ,  Z  ) .  Then,  (i)  T(«,Z)  -  «[W(Z)]  is  NPDF 

wrt  f l*  for  each  (measurable)  function  ♦(•);  and 

(ii)  G,  S(Z))  ■  t|i* [G,  S(Z)],  when  Z  is  governed  by  G,  is 

%  ^  % 

PDF  wrt  Jl*  for  each  (measurable)  function  <(»*(•)• 


S.  The  K-S  (Kolmogorov-Sairnov) .  Lilliefors  and  Srinivasan  Statistics 

(A)  Kolmogorov's  Original  Statistic 

Kolmogorov  (1933)  introduced  the  K-S  statistic  Dn(F0)  * 
sup|F  (z)  -  Fn(z)|,  for  continuous  cdfs  Fn(*)  and  empirical  cdfs 


F  (•).  F  (z)  *  -  J  e  (*  -  X.)  -  -  l  e(z  -  X(j)), 

n  n  n  j-1  3  n  3-1 

u  ^  0;  and  -  0  if  u  <  0. 


and 


e(u)  -  1  if 


Definition  S.l.  If  Xj,  ...»  Xn  are  i.i.d.  FQ,  continuous,  then 
Dn(F0)  *  K-S(n) . 

In  order  to  apply  the  K-S  statistic  directly,  one  must  know 
Fq(-)  completely.  However,  in  many  signal  detection  problems,  FQ(  • ) 
is  known  only  up  to  a  nuisance  parameter,  or,  equivalently  is  known  only 
to  be  a  member  of  a  specific  (parametric)  family .  Lilliefors  (1967,  1969), 
Srinivasan  (1970),  and  Choi  (1980)  introduced  modified  versions  of  the 
K-S  statistic  for  such  situations. 


(B)  Lilliefors 


Statistics 


Let  fl*  ■  (F(8 ;  • ) :  8  e  <0}  be  a  famil y  of  cdfs  admitting  a 

MLE,  8  ■  S(Xj,  ....  Xn),  for  8. 

Definition  S.2.  (1)  P(*)  is  the  cdf  satisfying  F(z)  ■  F(§;z) 

for  all  z;  and 

(ii)  6n  -  sup|Fn(z)  -  F(z)|. 
z 

Li 1 lief ora  (1967,  1969)  calculated  Monte  Carlo  tables  for  6  in 

n 

the  normal  and  exponential  cases,  while  Y.  Oioi  (1980)  has  given  such  a 
table  in  the  uniform  case. 

Srinivasan  (1970)  replaces  F(>)  with  the  Rao-Blackwell  estimate 
of  F(8; •)  in  the  Lilliefors  statistic. 

(C)  Srinivasan- type  Statistics 

Consider  a  family  (1*  »  (F(8;*):  8  e  ©}  of  cdfs  admitting 

a  M-S-S,  S(Z)  for  8. 

Definition  S.3.  (i)  F(z)  -  <  z|S(£)};  and 

(ii)  Dn  -  sup|Pn(z)  -  F(z)l. 
z 

Srinivasan  (1970)  computed  critical  values  of  Dn  by  (Monte  Carlo) 
sianilation  for  the  exponential  and  normal  families.  Some  of  his  numerical 
results  however  were  in  error,  as  was  pointed  out  by  Schafer,  Finkelstein 
and  Collins  (1972). 

Remark.  These  three  statistics  are  in  many  cases,  asymptotically 
equivalent.  Calculations  of  Fn(z)  for  the  case  when  A  is  unknown, 
s  ■  sQ  known  and  the  case  when  both  A  and  s  are  unknown  are  presented 
in  Appendix  A.  These  statistics  are  summarized  in  Table  5.1. 


Smirnov,  Lilliefors,  Srinivasan  Statistics 


TABLE  5.1.  (Continued) 


On*  can  now  treat  the  detection  p rob lee* . 


6.  One-Saaple  Detection  Procedure* . 

In  this  section,  one  will  derive  detectors  for  deciding  between  PN 

and  the  alternative  N  ♦  S  (noise  plus  signal).  The  data  set  will  be 

denoted  by  Z  •  (X.,  X,,  ....  X  )  and  a  is  the  PFA.  There  are 
i  »  n 

altogether  8  problems  considered.  Hie  first  4  deal  with  the 
cases  when  at  least  one  of  A  and  s  is  known  and  the  last  4  con¬ 
cern  cases  in  which  both  A  and  s  are  unknown.  Niawrical  examples 
for  each  case  are  provided  in  the  Appendix  B,  where  they  are  nurtured 
the  sane  as  the  cases  they  illustrate.  The  results  in  this  section  are 
sianarized  in  Table  6.1  at  the  end  of  this  section. 

Case  6.1.  (A  unknown,  s  ■  sQ  known,  l(A,s0)  t  Q(PRP)) 

PN:  A  ■  Aq  vs.  N  ♦  S:  A  +  Aq 

The  ainiaun  PDF  procedure  is  to  decide  N  ♦  S  if  and  only  if 

-1/ns- 

X(l)  <  Ajj  or  X(l)  >  bAg  where  b  ■  a  .  equivalently,  one  can 


use  the  statistic 


which  has  a  distribution  under  PN.  Note  that  the  procedure  is  based 

solely  on  the  M-S-S  for  A,  which  is  A  ■  X(l).  The  one-sided 
detection  procedure  for  PN:  A  _<  Ag  vs.  N  ♦  S:  A  >  AQ  My  be 
formulated  similarly. 

Case  6.2.  (A  unknown,  s  •  sQ  known;  class-fit  problem) 

PN:  l(A,s0)  e  Q(PRP)  vs.  N  »  S:  L(A,sQ)  i  Q(PRP) . 

One  employs  the  K-S  statistic  with  site  n  -  1  (K-S(n  -  1)) 
through  anyone  of  the  four  versions  of  N-S-N,  (V^>,  (U^>,  {N^},  {Y^> 

given  in  Theorem  3.5  and.  Definition  3.3.  Expllcitely,  they  are 


/n  .  n-1  -sAt 

■  sup  j  —r  l  e  (t  -  V.)  -  (1  -  e  0  )| 
*  z>0  j-1  3 

D<2)  *  ,up  I  JTT  c  <*  *  ui>  ”  O  -  *  *°)l 

n  t>l  n’1  j-l  3 

«£3)  "  ,up  I  iTT  1  e  (t  -  N.)  -  t| 

n  0<t<l  "  1  j-l  3 


i  n‘1 

»*»  I  I 

0<x<l  n  1  j-l 


e  C*  -  Yj)  -  t| 


where  *  JC-S(n  -  1)  for  i  •  1,  2,  3,  4.  Therefore  one  decides 

n 

N  ♦  S  if  and  only  if  >  d'  where  d'  is  an  appropriate  critical 

value  from  the  K-S(n  •  1)  distribution. 

Alternatively,  one  may  use  the  Srinivasan-tvre  statistic  (see  Table 


as  follows. 


- 


i 

i 


D_ 


(2-y)  *up 
"  1<  u  < 


.  ^  "i!  * tu  •  v  •  “  • u  ’0)| 


Not*  that  D  ■ 

n  n  n 


j 


The  Lilliefors-type  statistic  is  asymptotically  equivalent  to  the 
Srinivasan-type  statistic  in  this  case. 

Casa  6.3.  (A  •  A^  known,  s  unknown,  L(Aq,s)  e  Q(PRP)) 


PH:  s  ■  sQ  vs.  H  ♦  3:  s  i  sQ 
n 


The  M-S-S  for  s  is  **n(  \  in  XP^]"1and  the  statistic  X(j)  is 

.  J-l 

T  »  2sa  1  In  %  xL  under  PH.  Thus  the  detection  procedure  is 
0  Jil  *0  "»• 

•*ciu  "*•  itmtmtru  t  >  «  T<  x^.o/jr 


Hots  that  T  depends  on  Z  only  through  s.  The  one-sided  procedure 

% 

tor  PH:  s  ^  sfl  vs.  H  ♦  S:  s  >  sfl  nay  be  formulated  similarly. 


Case  6.4.  (A  »  Ap  known,  s  unknown;  class-fit  oroblea) 

PH:  L(Aq,s)  c  Q(PRP)  vs.  H  ♦  S:  t(A0,s)  i  Q(PM»). 

Again  we  should  use  the  K-S  statistic  with  the  M-S-N  which 

T1  T2  Tn-1 

according  to  Table  2.1  is  {y-  ,  y  ,  ...»  y  —)  +  U-0-S(n  -  1). 

n  n  n 

l  n"1 

Thus  let  D_  •  »vf>  |-f?-  I  e  ( t  -  yO  -  x  |  and  one  decides 
*  0<*<1  "  1  j«l  *n 


I 


.  24  . 


8 


I 


N  ♦  S  if  and  only  if  0.  >  d*  tdiere  d'  is  the  appropriate  critical 

n 

value  fro*  the  K«S(n  -  1)  table. 

Case  6.5.  (A,  s  both  unknown,  l(A,s)  c  Q(PRP);  Goodness -of- fit  test) 
PN:  l(A,s)  •  *-CA0.sq)  vs.  N  ♦  S:  l(A,s)  i  IfA^) 

Fro*  Theore*  3.3  (vi)  and  wider  PN  situation, 

h^br1 0 . ,w'1  *' ' 

Hie  statistic  is  (as  in  Case  6.3). 

T  •  -2  J  0  %  x^,,  and  the  decision 

rule  is  decide  N  ♦  S  if  and  only  if  T  >  Xfcn,,-*/!)  •*“  T  «  ^20.0/2)* 

Case  6.6.  (A  and  s  unknown,  l(A,s)  e  Q(PRP)) 

PN:  s  £  sQ  vs.  N  •  3:  s  >  sQ 

The  detection  statistic  here  should  only  involve  the  M-S-S 
S(2)  •  s-n(  l  In  IWI]'1  for  s. 

%  j«|  MU 

2*0 

The  decision  rule  is:  Decide  N  ♦  S  iff  -7-  <  c*,  where  c’ 

s 

2 

is  the  (100cs)th  percentile  of  the  X(2n.2)‘**#trtbut*0,,‘ 

One  notes  that  if  *  In  X^ ,  then  the  Y's  are  i.i.d.  T-exp 
(In  A,s).  Park  and  11.  Choi  (1978)  derive  the  nininun  PPD  one-sided 


1 


procedures  for  Che  shape  parameter  s.  It  will  be  shown  in  Appendix 
A. 3,  that  the  Park-Choi  decision  rule  is  equivalent  to  the  one  given 
above. 

For  the  two-sided  detection  p rob lea:  PN:  s  »  sQ  vs.  N  *  S: 

s  t  sQ,  the  decision  rule  follows  froa  that  above. 

Case  6.7.  (A,  s  unknown;  l(A,s)  e  H(PRP)) 

PN:  A^A0  vs.  H  *  S:  A  >  AQ 

One  should  employ  the  H-S-S,  (A,  s),  and  the  statistic 

T(Z,  A)  •  rJjUffJQ)-!.,1 S-Al  -x.  F(2,  2n  -  2). 

V  f  In  1.  •  n  In  X(I) 

j-1  5 

The  decision  rule  is:  Decide  N  ♦  S  iff  X(l)  >  AQ  and  T(Z,  AQ) 
where  f  is  the  100(1  -  o)th  percentile  of  the  F(2,  2n  -  2)- 

distri button. 

Again,  aa  in  Cue  6.6,  Park  and  0.  Choi  (1978)  give  a  ainiaun  PDF 
procedure  for  the  truncated  exponential  cue,  which  is  equivalent  to 
the  decision  rule  above. 

Cue  6.8.  (A,  s  uiknown;  clus-fit  problem) 

PN:  l(A,s)  G  O(PRP)  vs.  N  »  S:  l(A,s)  fc  O(PRP) . 

One  cm  ac  the  Kolaogorov-Saimov,  Lilliefors  or  Srinivuan-type 


« 
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statistics.  From  Theorem  3.4,  one  finds 


-  sup  ] 
n  0<z<l 


D<2>  - 


i  n"2  Vi 

“2  i  e  u  -  -y — )  - 

n  j-l  n-l 


1  n;2  .  * 


n  K  1  n-2  t 

n  0<z<l  j-l 


-1  n-l 


then  K-S(n  -  2)  for  i  ■  1,  2.  The  decision  rule  is  decide 

n 

N  ♦  S  if  and  only  if  D^>  d'  where  d'  is  the  critical  value  from 

n 

the  K-S(n  -  2)  table. 

One  now  considers  the  Srinivasan  statistic.  From  Table  5.1,  one  has 

Dn  "  SUP  lFn(2)  ' 

•  SUP  n/s’!n  I  €  (XTTT  *  Ym)‘1^nT')l1'  V  ln  a1"’11 
a  <  z  <  ae  j-l 

where  A  -  X(l)  -  a,  s  -  s'.  Let  v  -  yyjy  >  v(j)  ■  •  Th«n 


f\j  i  I  "  * 

d  -  sup  /,!—♦—  [  y 

"  i<7<  ,n/I  n  »  jii 


i  .  i  (  j  £  (v  -  V())))-l  ♦  -  jj  In  *]"■'' 

n  n  ..  n  n 


*  .  »“P  „/,•  Is  IT  E  »"  v  -  *"  *  'T1'"  •  n  »"  v]”'11 


n/s'  'n  1  L 
1  <  v  <  en/5  n  j-l 


Let  u  -  ln  v,  U(j)  «  ln  V(j),  then  since  s  -  ■  -t5 - -  - 

i  m  V4 
j-l  J 


S'  ..ntt-l. 


has  D  -  sup  |{-r  [  £  (u  •  U(j))-  {1  -  [1  -  J  u] 

n  n  0  <  u  <  n/s'  n_1  j-l  " 


(~)  sup  _  l~r  n'l  e  (u  -  U(j))-  (1  -  [1 - “j"  S| 

0  <  u  <  (n-l)u  n_1  j-1  (n-l)u 


•  SUP  n-1  ICT  A  G  <u  *  U<3»  *  {1  *  »  " 

0  <  u  <  l  U.  3-1  (n_1)u 

j-1  3 

The  critical  values  of  this  statistic  cannot  be  obtained  from  the 
known  existing  tables.  They  may  be  obtained  by  the  Monte  Carlo  simulation 
method  for  various  sample  sizes  and  PFA 

These  statistics  are  summarized  in  Table  6.1. 


7.  Two -Samples  Detection  Procedures 


The  data  set  here  is  Z  «  (X,Y)  ■  (X.,  ....  X  ,  Y. ,  ....  Y  )  - 

%  %  1  mi  II 

(Z, ,  Z_,  ....  Z*.)  where  N  »  m  ♦  n.  Here  X,  Y  are  two  independent 
random  samples  from  Pa(Aj,  s^)  and  Pa(A2,  s2)  respectively.  The 
letter  a  will  denote  the  PFA.  As  in  the  one-sample  case,  numerical 
examples  are  given  for  each  case  in  Appendix  B  with  the  same  order 
and  numberings  as  they  are  presented  here. 

Case  7.1.  (Aj,  A£  known,  Sj,  *2  unknown) 


PN:  Sj  -  $2 


vs.  N  ♦  S:  Sj  i  s2 
m  X, 


n  Y. 


The  M-S-S  for  (s.,  s,)  is  (s. ,  s,)  «(m(  \  ln(-j^-)l  l,  n(  l  lnCr3-)]  *) 

12  12  j-i  Al  j-1  a2 


TABLE  6.1.  (Continued) 


and  under  PN,  T*  = 


^  F(2m,  2n) .  Therefore, 


decide  N  +  S  if  and  only  if  T  >  f*  or  T  <  f"  where  f*  and  f" 
are  the  appropriate  percentiles  of  the  F(2m,  2n) -distribution. 

An  important  point  arises  in  Case  7.1  above.  The  detection 
statistic  T*  involved  the  M-S-S  (s^,  s^)  for  (s^,  S2) ,  the 
unknown  parameter  pair.  However,  the  remark  following  Example  4.1, 
and  the  cases  of  Section  6,  suggest  that  for  Case  7.1,  one  should  employ 
the  M-S-N,  since  the  particular  values  of  s^,  and  s2  are  not 
pertinent  here.  This,  indeed,  is  the  case,  as  will  be  seen  from  the 
derivation  to  follow. 

One  can  directly  verify  (from  Lemma  3.1) 


Z. 


u  •  u  • 

Theorem  7.1.  Let  W.  =  In  for  1  <  j  <  m;  and  *  In  -r^-  for 

-  j  Aj  a2 

r  N  -1 

m  +  1  <  j  <  N;  and  V  =  [£  W.]  [£  W.]  for  1  £  r  N  -  1.  Then 

r  1  3  1  1 


(i)  (Wj,  ....  WN)  are  i.i.d.  Exp(s). 

(ii)  (V^,  •  •.,  VN  j)  is  the  M-S-N;  and  is  a*  U-0-S(N  -  1).  Further, 


(iii)  T* 


fHi  (JB _ 1 

'•nr  4-V  ] ' 
m 


This  means  that  T*  is  a  function  of  the  data,  Z,  solely  via  the 

% 

M-S-N  of  the  combined  sample.  Hence,  T*  is  both  a  function  of  the 
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a 


M-S-S's  for  the  individual  samples  and  the  M-S-N  of  the  combined 
samples. 

The  one-sided  detection  problems  may  be  handled  analogously. 


Case  7.2.  (A^,  A2  unknown,  s^,  s2  known) 


PN:  Ax  =  A2 


vs.  N  +  S:  Aj  i  A2 


A  A 


In  this  case,  the  M-S-S  for  (Aj,  A2)  is  (Aj,  A2)  *  (X(l) ,  Y(l)). 
Since  ln(— j-^-)  ExpCmSj),  ln(^^-)  ^  Exp(ns2),  one  has  by  Lemma  3.6 

the  following. 

Lemma  7.2.  Let  T  »  *  Then  unc^er  PN»  the  distribution 

of  T  is  given  by 


r 


ms. 


H(t) 


1  - 


(mSj  +  ns2) 


ns. 


(ms.  +  ns_) 


-ns  t 

e  t  >  0 


ms^t 


t  <  0 


The  proof  of  this  lemma  is  given  in  Appendix  A.  Thus  the  decision  rule 
is:  Decide  N  +  S  if  and  only  if  T  >  or  T  <  C2  where  Cj,  C2 
are  determined  by  H(Cj)  *  1  -  ^  »  a°d  H(C2)  =  y  . 

In  the  spe.cial  case  where  =  s2  =  sQ  is  known,  the  decision 
rule  reduces  to:  Decide  N  +  S  iff  Y(7)"  <  bj  or  >  b2  where 

bl  -  (f)]"80  andb2-  llfy]™0  • 


-  32  - 


One  further  notes  that  X-s  ^  £1(0,6)  with  9  =  A  s.  For  this 

uniform  case,  Y.  Choi,  Bell,  Ahmad  and  Park  (1982)  present  a  detection 
procedure  which  coincides  with  the  special  case  above. 

Case  7.3.  (Aj  *  A2  =  A,  unknown;  s^,  s2  unknown) 

PN:  si  3  s2  vs.  N  S:  s^  i  s2 

One  first  attempts  to  base  the  decision  rule  on  the  theorem  below. 

Theorem  7.3.  Let  X, ,  . . . ,  X  ,  Y, ,  . . . ,  Y_  be  i.i.d.  P(A,s). 

. . - —  l  n_i  n  l  n  n_j 

Then  (i)  X(l) ,  Y(l) ,  ^  In  and  n2  =  j  In 

are  independent;  with 

(ii)  2ns  |  in  X(l)  -  In  Y(l)  |  ^  X*  ;  and 

d  2 

(iii)  Zsrij  =*  2sn2  ^  X2n_2 

The  decision  rule  for  this  case,  with  m  =  n,  would  then  be: 

Decide  N  ♦  S  iff 

.  2ns |n,  -  n2|  In,  -  n2| 

q  3  - i - i -  »  — i - = -  >  c* 

2ns | lnX(l)-lnY(l) |  | In  X(l)-lnY(l)| 

However,  the  cdf  of  Q  is  not  known,  and,  hence  orie  seeks  other  approaches 
Beg  (1980)  derives  a  uniformly  minimum  PFD  procedure  for  the  case: 

PN:  Sj  <  $2  vs.  N  +  S:  Sj  >  s2 

in  the  truncated-exponential  case,  which  applies  to  the  Pareto  oroblem 


at  hand  even  when  m  +  n. 


ill  11 

Let  n.  and  n,  be  as  above;  n  !  I  In  X.  +  £  InY.-NW, 

1  3  j*l  3 

where  N  =  n  +  m  and  W  *  min(ln  X(l),  In  Y(l)}.  Beg  proves 


Lemma  7.4.  The  conditional  density  of  t^,  given  W  and  n,  is 
(m-m-2)!n™~2  (n  -  nJn_1 

h(nJw,n)  * - -  for  o  <  n.  <  n. 

1  (m-2) ! (n-1) !  n 1 


The  decision  rule  based  on  this  lemma  becomes: 

n 

h)dn1 


<  c  *  c(w,r|)  where 


rn 

,  j  h(Hj  |w. 


'1 


a. 


Decide  N  ♦  S  iff 


If  one  performs  the  actual  integration,  it  is  easily  seen  that  c 
satisfies  the  relation  c  »  c'z  where  Ic,  (m  -  1,  n)  *  1  -  a, 

-  y)n"1  * 

is  the  incomplete  Beta- function.  A  table  for  the  incomplete  Beta 
function  has  been  tabulated  by  K.  Pearson  (1934) . 


Case  7.4.  (Aj,  A2  unknown, 
PN:  Aj  <  \2 

Let  hj,  n2,  W,  n  be  as  defined 


=  s2  *  s  unknown) 
vs .  N  +  S:  Aj  >  A2 
in  Case  7.3.  Let  N  *  m  +  n. 


h(x*|w,n) 


if  x*  *  m 


SBffiil  [1  .  lG£2l.]  if  »  <  X*  <  w  ♦  3  , 

Nn  R  ra 


and  define  the  number  c  ■  c(w,  n)  by 


h(x*|w,n)dx*  ■  a. 


From  the  result  of  Beg  (1980) ,  the  decision  rule  is  decide  N  ♦  S 
if  and  only  if  In  X(l)  >  c.  The  number  c  may  be  found  by  performing 
the  actual  integration  in. which  case  one  gets 

1 

n  ,  ^  fT2 

c  a  +  w  where  c'  *  1  -  (^) 
n  n 


Remark :  Suppose  m  *  n  *  y  ,  then  Lemma  3.6  (vi)  may  be  used 
to  test  PN:  Sj  ■  s2  vs.  N  »  S:  Sj  f  s2  in  Case  7.3  and  to  test 
PN:  Aj  ■  A2  vs.  N  +  S;  Aj  f  A2  in  Case  7.4.  In  both  cases,  under 
PN,  the  statistic 

T  .  (N-2)»|ln  X(l)  -  in  Yjlli  %  p(2  2N  .  4). 

4(h1  ♦  n2) 


The  decision  rule  is  decide  N  +  S  if  and  only  if  T  >  f'  where  f' 
is  the  appropriate  percentile  of  the  F(2,  2N  -  4) -distribution. 


Case  7.5.  (Aj,  Aj,  Sj,  s2  all  unknown) 


PN:  -  s2 


vs.  N  +  S:  Sj  *  s2 


Let  n 


then  from  Lemma  3.6  (iii),  under  PN, 


(n  -  1 )Hj 
(®  -  i)n. 


^  F(2(m  -  1),  2(n  -  1)). 


The  decision  rule  is  decide  N  ♦  S  if  and  only  if  T  >  f*  or  T  <  f" 
where  f'  and  f"  are  appropriate  percentiles  of  the  F(2m  -  2,  2n  -  2)- 
distribution. 


Case  7.6.  (Aj,  A2>  si>  s2  al*  unknown) 

PN:  Aj  ■  A2  vs.  N  *  S:  Aj  i  A., 

This  detection  procedure  does  not  appear  to  have  an  elementary 
solution.  One  may  try  to  apply  the  likelihood-ratio  test.  One  has  when  not 


assuming  Aj  »  A2> 


\  • x(i)’  s~i  * 


A2  -  Y(l),  s 2  *  n[  In  . 


Under  PN,  one  has 


A*-  min  (X  (1),  Y(l) }  -  Z(l)  and 
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if  Z(l)  -  X(l)  <  Y(l) 


if  Z(l)  »  Y(l)  <  X(l) 


The  distribution  of  L^/L Q  is  complex  and  the  above  expression  resembles 
in  structure  that  of  the  Behrens -Fisher  problem. 

Case  7.7.  (Aj,  A2 ,  Sj,  s2  all  unknown) 

PN:  (Aj ,  Sj)  ■  (A2,  s2)  vs.  N  ♦  S:  (Aj,  Sj)  i  (A2>  s2) 

The  decision  rule  consists  of  2  steps,  first  deciding  whether 
Sj  *  s2  and  then  if  one  decides  "  *2  one  t**®8  to  decide  whether 
A^  *  A2<  The  procedure  is  a  combination  of  Cases  7.5  and  7.4. 

ut  ni  ■  jl, 10  •  nz  ■  J, ln  •  then 

(n-l)n. 

under  PN,  T  -  j-  .x~  *  F(2(m  -  1),  2(n  -  1)). 

v®“  *■ )  *1 2 

The  decision  rule  is: 

(i)  Decide  N  ♦  S  if 

T  >  f(2(m-l) ,2(n-l)l-a/4)  5  C1  or  T  <  f(2(m-l) ,2(n-l) ,a/4)  5  C2 


(ii)  If  C2 
In  X(l)  >  C3 


<_  T  <  C. ,  then  decide  N  ♦  S  if  and  only  if 

1  nc* 

and  In  Y(l)  >  C.  where  C,  ■  — -  ♦  w  C  • 

4  3  n  4 


w  «  «in  {In  X(l),  In  Y(l)}, 


C3 


>  .  &-2 


C4 


■  n 

n  -  l  In  X.  ♦  l  In  Y.  -  I 
j»l  5  j»l  J 


Na)N- 


N  ■  ■  ♦  n. 


Table  7.1  smaarizes  these  procedures. 


TABLE  7.1.  Two-Saaples  Signal  Detection  Procedures 


TABLE  7.1.  (Continued) 


< 
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I 
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APPENDIX  A:  Proofs 


1  .  (Hli)  (  *Jil)  2  >  x(i) 


n  ' '  z 


A.l:  A  unknown,  s  =  sQ  known,  then  Fn(z)  = 


^  1  X(l). 


1 


Proof.  Fn(z)  =  E{I{X  <  2}(SCZ)}  =  P{Xj  <  z|x(l)} 


'1  - 


=  £  I  PCX.  <  z|X(l)> 
3*1  J 


I  P{X(j+l)  <  z|X(l)}  z  >  X(l) 
3*1 


n-1 


n-1 


“  n  *  n  %  P  —  xflj  from  Theorem  3.5  (ii) 


3*1 


Now  let  U(j)  =■ 


Jlhll 

X(l) 


,  then  (U(l),  13(2) ,  ...,  U(n-l)}  *  Pa(l,  sQ)- 


0-S(n  -  1). 

Therefore  F  (z)  =  —  ■  'n"1' r  1 
n  " 


xfn 


)] 


1  .  .n-1. 


-j* 


1  -  PtV  w’ 


i  -  (— i  (itii)  . 


n  '  '  z 
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A. 2.  A,  s  both  unknown,  then  Fn(z) 


z  <  a 


z  =  a 


(— )  [1-  -  In  -]  a  <  z  <  ae 
v  n  *  1  n  aJ 


^  ^  n/s 

<  z  <  ae 


where  A  =  a,  s  *  s. 


Proof.  Fn(z) 


E{l{Xl  <  z}lS^}  =  <  z|CA,  s)} 


~  l  P{X.  <  z|(A,  s)} 

j-1  3  ~ 


~  ~  l  P(X(j+l)  <  z|(A,  s)>. 

"  "  j  =  l 


z  >  a 


Let  U(j)  *  ,  V(j)  *  In  U(j),  then  by  Theorem  3.5, 


{U(l),  U( 2),  ....  (J(n  -  1)}  'v  Pa(l,  s)-0-S(n  -  1)  and 


{V(l),  V(2) ,  ....  V(n  -  1)}  -v  Exp(s)-0-S(n  -  1) . 


Hence  F  (z) 
n 


-  -  ♦  -  nf  p{ua)  <  fis} 

n  n  .fc,  —  a1 

3*1 

=  -  +  -  V  p{V(j)  <  In  |-!s> 

n  n  .L.  J  —  & 

3*1 

,  .  n-1  V.  _  _ 

.I  +  I  l  P{~J— <1  Inf) 

n  n  jil  (n-l)V 


since  s  *  s  * 


n-1 

l  Vj 

j-1  3 


(n-l)V 


- 


I 

I 

! 

! 

! 
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v  V  ♦  V 

By  Lemma  3.1,  { - —  , 


(n-l)V  (n-l)V 


,  . ^  U-0-S(n  -  1).  To 


continue. 


F  (z)  =  -  +  (— )  P{ - - —  <  -  In  -} 

n  n  n  (n-l)V  "  n  a 


=  1  -  (~)  P  { - - -  >  J  In  i> 

n  (n-l)V  n  a 


"  1  -  C^Hl  -|ln  f]“'\ 


In  the  above  derivation,  since  <  (n  -  1)V 

may  restrict  z  to  In  —  <  —  or  z  <  aen^s. 

a  —  s  — 


n 

s 


for  every 


one 


A. 3.  One  verifies  here  that  the  u.m.p.  procedure  found  by  Park  and 
U.  Choi  (1978)  is  equivalent  to  that  in  Cases  6.6  and  6.7.  Park  and 
Choi  considered  the  p.d.f.  of  the  truncated  exponential. 


f(x,A,v)  »  I(v,»)  ^ 

where  A,v  are  unknown.  Let  Y^,  Yj,  .. 
a  truncated  exponential  distribution.  Let 
Y(l)  =  y. 


0  <  V,  A  <  ® 


. ,  Y^  be  random  sample  from 
n 

S  =  £  Y. ,  then  given 

i=l  1 


f(s|y) 


An_1 (s -ny)n-2  e-X(s-ny) 

r fcT) 


1  (ny ,°°)  ^ 


(1) 


while  given  S  ■  s 


4' 


I(v,s/n)  (y)- 


Theorem  A.  (Park  and  U.  Choi).  For  testing  hypothesis  Hy:  v  <  v 
against  Ky:  v  >  vQ,  the  u.m.p.  unbiased  test  is  given  by 


<KY(1))  = 


1  if  Y(l)  >  C(s) 


0  if  Y(l)  <  C(s) 


where  C(s)  is  uniquely  determined  by  P(Y(1)  >_  C(s)|S  =  s,  v  =  vQ}  =  a 

Theorem  B.  (Paik  and  U.  Choi).  For  testing  hypothesis  H^:  X  £ 
against,  X  >  XQ,  the  u.m.p.  unbiased  test  is  given  by 

(  1  if  S  <  C(Y(1)) 

<KS)  -  \ 

I  0  if  S  >  C(Y(1)) 

where  C(Y(1))  is  uniquely  determined  by  P{S  C(Y(1))|Y(1)  *  y,  X  *  XQ)  =  a 
From  Theorem  A  and  (2),  one  concludes  C  *  C(S)  is  determined  by 


F  (s  -  ny)n-2  dy 
(s-nvj"-1  'C 


(n  -  1)  f1  (1  -  u)n"2  du 
JC* 


n(Y(l)  -  v  ) 

where  U  =  - 5 - — - -  ^  B(l,  n  -  1).  But 

s  -  nvQ 


n(Y(l)  -  v  )  T, 

U  ■  S-.Y(l)^Y»)-nv0  *  TTrY  wher,  V  ■  (n  -  l)T. 


8 


T  is  the  statistic  given  in  Section  6,  Case  6.7.  Since  the  function 

t1 

f(t')  =  -+- p-  is  strictly  increasing  for  t'  >  0,  the  two  tests  are 

equivalent.  It  is  also  easily  verified  that  if  U  ^  B(l,  n  -  1) , 

T* 

then  — - j-  ^  F(2,  2(n  -  1)).  The  proof  of  the  equivalence  between 

Theorem  B  and  the  test  procedure  in  Case  6.6  follows  from  (1)  in 
the  same  manner  and  will  be  omitted. 


A. 4.  (Proof  of  Lemna  7.1).  From  the  hypothesis,  let  U  *  lnY(l), 


V  =  In  X(!)  and  T  *  U  -  V,  then 


P{T  <  t}  =  P{U.  -  V  <  t) 


-ms  v 

I  P{U  <v+t}(ms.)e  dv 

Jo  - 

E-ns~(v+t)  -ms..v 

[1  -  e  ]e  dv. 

ax(-t,0) 


Considering  separately  the  case  t  >  0  and  t  <  0,  the  above  integral 
can  be  evaluated  and  equals  to  H(t)  given  in  Lemma  7.2. 


APPENDIX  B:  Numerical  Examples 

The  numerical  examples  correspond  to  the  various  techniques 
developed  in  the  main  body  of  the  report,  and  have  corresponding  numbers. 
The  data  is  divided  into  three  sets. 

Table  1.  This  consists  of  simulated  interarrival  times  from 
Pa(l;l);  Pa(2;l);  Pa(3;l);  and  Pa(5;5)  distributions. 

Table  2.  This  consists  of  simulated  interarrival  times  from 
Pa(l;2);  Pa(l;5);  Pa(2;2);  and  Pa(5;2)  distributions. 

Table  3.  This  consists  of  real  data  related  to  24  complete 


heart-beat  cycles.  [Special  attention  will  be  paid  to  the  waiting 
times  for  the  "R-peaks".  See  sketch  below  the  table.] 


Pareto  Data  (Interarrivals) 
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«4,0«n<vom«i*)oao9'9>O'»4^*i«<'i**i/'N<ro«^-«’<otf\^fw«4r»4>o<*t<nO‘<OiftO'6'r«>oi>*e*4’<r4O‘ 
OONN*n«'0N<,>0'fi®®(<N^I««N>J»0lft00l,l<,<O»'<fN0,,<0^#^N^tACWMOO'#O®< 
^  »#»«••••••••••••••••••••*•••••••••*••••••••••••••• 

to  «r»<r»  <0  **  ifM/MA  « r\j  »/N  ^  ia  kma  or*  i*' 'Ox1' '‘*“,'Oirv  if*  u'if\r-»r‘»‘>o 

• '  ^  0*  prf  H  H 

l/> 
w 

a. 

>4><«4>4«4i4WH«4HNiMNrwiwNNiMiyiNnmmnctnmnnn«ee^4r*eeeeir 

N0^0N<MM*4N«N4>oNHNnseOHei4o«on>40KmN0eo>>«iee>4«e>N»e#«'e0o«4ie 
«<4OiANK«H»«»«H^w»tf\Ni4oe«9N0NNK»omom«iivn«om«o#e>4#A<««eo 
0NN0Ka^NlA9Hnn0NmN0n99«NNNO9N^^>4O««^O«4NU>nHn>«OO<^N#OA 
/->  ••*•••••••••••••••••••■•••■••••••••••••••••••••• 

1-1  0«n0w«n*4ie>4»«4>4en09iA»AnAe»«neoeNm«nonNnn«>4nnoin«aNn» 

«  N  *4«4p4fM  «  w*  m4  m4  m4  m*  *1  «  N  N 


^^•(x^HHHHNNNNNiMNNNNnnmninnAmnmeeeeeeeeeetfN 


0004«ivOoo>400m<eo«o<t0iAnN^#^^oN0>4O0Oieoose'^0«<e<MiMMmH0ON 

#«4»0N»^a#ONON«4N»«No«NO»ee»o«N#irANmninNiAO«o«««o>4mf)»HA 

^e#«a^in<v««Heirecsoo^N«iNMNeN^ineNnoM<'ir»>4«NN«NOO(Moeo>4« 

•••••••••••••••••••••••••••••••••••••••••••••••»•• 

nN«NinNiM0#(n«<rnNennmNnNNmN«NinO'n#iM«0mKMANnmienNNO«Nn>««« 

-CM  •*  « 


*4Mneie«N«»e^Nn'fir<0^«»OHNneia«Na»o^Nn«i<\«N«»o>4N<'<«««N«»o 

»4>4M«<w>4(4»4iH«4ivMNNNNNNNNnAAfiru^n«nmme#e#eeeeee^ 


AIMNrte  »  CO^lAOON^NeN  A»A*«VlA^I»»®IA»«4»lf'»V 

NONmN«4N#«ieNNH««roo<os«HNe>4i4neeNH«io^«4O0N"«^N^NAeNO 

•4«4HD>N»nm«4NN>4Hen>4*(>4«4<rf«4«N««M<Mmn<4Hn^p4>4NIV>4Nf4«MN«4NNe^lMM 
**  *4  «4  tft 


TABLE  2  Pareto  Data  (Interarrivals) 
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IN 
•  « 
m 


a 

a. 


iwn<»o‘f*0‘c4iv(v*0‘iv©»-4©rt^«'4^o*'*©«*«i4*o©©*»l''»*»vivfvo<<"i<rM.«*,fn©.r«ft<r5i©©«#*iv©© 

»#*##•#••••••••••••••••••••••••••••••••••••  •  •  •  •  •  •  • 

WN>49He4^>i>*4>4'A^4^^4eo«4NS4H«44MenN»e>4»ie(Meoo>^o«ie<e4ie 

N|4«4«4  *4  <V*4  N  H  H  W  9“  ^ 


H»4>4^*<H<4i4«4HNiMNNNNMNNNnn^n(nnninn«'e#«ee#e«#« 


«n~<e»tfit©~4**04ro»©tfN©^««©Nm»4<n<vfvm4B«*4«'«nK<MM>‘©0‘<c^©«~4tfM*>aMfwv'«>*»'Pi‘© 
oNnn«>4««>«SHMei0ienNH#mnefg»>4NNan«<4B0eo«««0ONn90Ni«iie<4e^e 
4«^M«<4^H4,0#4O>4NNNON4inn«e«4NN>4ipniMOM«NooomieNO<4mNO‘Ni4e'e 
/— >  •••••■•••••••■•■•■•••••••••••••••••••••••••••••••• 

N  nMVN0><VNNNN#«mNieeNNNNNnNAiNnNNieS)VIMniNNNNNIMmNIMNmNNNNinm 

•  4* 

CM 


.«■»  *\»  ••*♦»«'*>  p-©<y-o*-«»M  t*» -e  «rs  «o^«ocr  ©-<  *»"'-*v'"O^«0‘© -*  im  «»,«■»/*  ■0^*0  <>o«*  cm '*i*-m,0^io^© 

,4*4  *4  *4  ~4  *4 -4  *4  .4  *4  <vrvrv«v<vfv<vfvfv<v<'im<n<n<n<*w*trw'> 


<0O(niennirv»m>4MM«9nNiea«nMP(nn^<r9«ie9«o«o»»>4«N<0>4.<>4>4iueiivM'e 

<AN4N^««»ONci»o»iA>4#e©«ie<4MnoninN^«M><#oieoNin^OHivMA«4neiAn 

©•4«o^o^0inrtNmno^o»4<4ONOMN«4>4ON)V<4>4om(nHOOoonMOOO«e'4n>4iA«<N 

M»4»4>4»<i  *<  »ml<  li'4ft|i4»4»<»<><t4«4»4>4»<»4><Hii  <■  <i  «»nCi  In4»  li  <•  «<  <■  li  <«  l*4<  li  Mp<  *  MH»<»4>4M>'< 


« 

a. 


CM 


•<Nn#ie«N09O«iNM«ir«K«»o«4Nn«ir<eN«»o<4Nn#ie«Ne»o^iMA#iA«Ne9O 

•4*4  •*  *4  »4 1>4  »4  «« «4>rfe4  MfsiNMN  iv  <w4Mfw<n<n<*t<n<*u''i  <*)<**<<*<**  «'4r*<r*4r*'r«*irt 

♦•©»i>©-»fn©©©ivp**iv-»©n©«n«©^"r|rr|»vi2©«'g-<r;^J{^~02SoJ2!lJS2*isi2S45Sln 


cd 

Ou 


*M*<^«<»M^*4^-«-4'Vfvivfv«vrsjpMMMMinr,w^<*'(^«nrvnc-iin4’4-*r  «■  >r 


TABLE  3  EKG-Data  (Heart  Beat  Cycles) 


-  52 


MAGNITUDE 


TIMES 


Cycle 


1  p 

Q 

R 

.0927 

.0708 

1.190 

.100S 

.0917 

1.203 

.1121 

.0982 

1.212 

.1238 

.0608 

1.297 

.0819 

.0206 

1.261 

.0769 

.0764 

1.284 

.0815 

.0507 

1.228 

.1263 

.1113 

1.166 

.1009 

.0817 

1.274 

.0839 

.0601 

1.239 

.1126 

.0907 

1.164 

.0939 

.0894 

1.189 

.0929 

.0648 

1.278 

.1010 

.0555 

1.213 

mm 

.0560 

.0374 

1.235 

.0832 

.0899 

1.294 

.1278 

.0939 

1.294 

.1044 

.0585 

1.255 

.1184 

.0702 

1.250 

.1142 

.1141 

1.266 

.1338 

.1140 

1.228 

.1331 

.1050 

1.230 

.0849 

.0584 

1.127 

.0916 

.0807 

1.156 

.1442 

.1323 

1.260 

T 


.4785 


.4344 


.4316 


.5652 


.5060 


.4757 


.  4  7C0 


.4240 


.4542 


.4663 


.4462 


.4554 


.4773 


.4742 


.4874 


.4263 


.4465 


.5001 


.4451 


.  4105 


.4052 


.4405 


.4265 


.402  3 


Cycle 

Length 


1014 


1014 


1040 


1041 


1056 


1078 


:  1073 


1076 


1089 


1044 


1012 


1028 


1038 


1062 


1063 


1083 


1063 


1018 


1013 


947 


891 


911 


901 


947 


Example  6.1:  A  unknown,  s  =  s^  ^nown»  £(A,Sg)  e  ft(PRP) 

PN:  A  ■  Ag  vs .  N  ♦  S :  A^AQ. 

Detector  Statistics :  (i)  X(l)  or 

(ii)  T  -  2ns  «  In  ^  . 

0  A0  2 

Decision  Rule:  Decide  N  ♦  S  iff 

(i)  X(l)  <Aq  or  X(l)  >  bAg,  b  -  V1/ns0 

(ii)  T  <  0  or  T  >  x|>(1.a). 

Decision  rule  (i) ,  the  MP  procedure,  will  be  illustrated  using 
the  following  PN  situation: 

PN:  A  «  2  vs.  N  ♦  S:  A  )l  2  (s  known,  a  ■  .01) 

Data  Sets :  (See  Tables  1  and  2) 


1. 

V 

•••»  V50 

i.i.d. 

Pa(l;l) 

2. 

V 

•••*  *50 

i.i.d. 

Pa(2;l) 

3. 

V 

X50 

i.i.d. 

Pa(3;l) 

4. 

V 

’ -  ’  Y50 

i.i.d. 

Pa(2;2) 

S. 

V 

....  z50 

i.i.d. 

Pa(l;5) 

The  Decision  Rules  are: 

1.  Decide  N  ♦  S  iff  V(i)  <  2  or  V(i)  >  2.193.  Since  V(l)  «  1.010 
one  decides  N  ♦  S. 


2.  Decide  N  ♦  S  iff  *f(l)  <  2  or  W(l)  >  2.193.  Since  N(l)  «  2.033 
one  decides  PN. 


3. 

Decide  N  + 

S  iff 

X(l) 

<  2 

or 

A 

X 

2.193. 

Since 

X(l) 

*  3.001 

one  decides 

N  +  S. 

4. 

Decide  N  + 

S  iff 

Y(l) 

<  2 

or 

Yd)  > 

2.094. 

Since 

Y(l) 

=  2.044 

one  decides 

PN. 

S. 

Decide  N  + 

S  iff 

Z(l) 

<  2 

or 

Z(l)  > 

2.037. 

Since 

Z(l) 

=  1.008 

one  decides 

N  +  S. 

Example  6.2 

A  unknown , 

s  « 

s0 

known 

PN:  L( A,sq}  e  fl(PRP)  vs.  N  ♦  S:  i.(A,s0)t  rt(PRP) 


Detector  Statistics:  Four  Kolmogorov- Smirnov  statistics 

(D^,  d„2^»  Dn^  *  Dn4^  811(1  8  Srinivasan-type  statistic  (Dn)  are 

available. 

Decision  Rules:  (1-4)  Decide  N  +  S  iff  Dn  >  Q  where 

d  ,  is  value  from  Kolmogorov-Smimov  table  and  D  ■  D^,  D^2^ , 
n-1  ,a  n  n  n 

°f3)»  Di4)-  (5)  decide  N  ♦  S  iff  D  >  (~M  .  „  . 

n  n  no  n^x^ot 


Data  Sets:  (see  Tables  1  and  2) 


p<? 


-  ss  - 


critical 

value  statistic 


statistic 


Data  Sets :  (see  Tables  1  and  2) 
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1.  Xj,  ....  X5Q  1.1. d.  Pa(l;2) 

2.  Yj,  ....  Yso  i.i.d.  Pa(5;5) 

Decision  Rule:  Decide  N  +  s  iff  T  >  140.2 
From  data  set  1:  T  =  91.017,  so  one  decides  PN. 

From  data  set  2:  T  =  48.737,  so  one  decides  N  + 

Example  6.4:  A  =  AQ  known,  s  unknown 

PN:  L(Aq,s)  e  Q(PRP)  vs.  N  ♦  S:  L(A0,s 

1  n_1 

Detector  Statistic:  D*  =  sup  I — r  T 
-  .  n  o  <  u  <  1  n-1  1 

r  '  X. 

where  Tr  *  I  IhCt1) 

j  =  l  0 

Decision  Rule:  Decide  N  +  S  iff  D*  >  d  , 
-  n  n-l,< 

is  appropriate  value  from  Kolmogorov- Smirnov  table 

Data  Sets:  (see  Tables  1  and  2) 

1.  Xj,  ....  X5Q  i.i.d.  Pa(l;2) 

2.  Y1#  Ysq  i.i.d.  Pa(5 ;5) 

Decide  N  ♦  S  iff  D|Q  >  .233 

Since  D£q(X)  *  0.101,  one  decides  PN. 

Since  D|q(Y)  *  0.0S3,  one  decides  PN. 


a* 


or  T  <  82.4 

S. 


£  «(PRP) 

T. 

e  (u  -  yi-)  -  u| 
n 


where  d  ,  _ 
n-1  ,a 
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( 

I 

! 


Example  6.5.  A  unknown , 
PN:  L( A,s)  =  l(Aq,Sq) 

Detector  Statistic  T  = 


s  unknown,  i(A,s)  e  O(PRP). 


vs.  N  *  S:  l(A,s)  +  L(A0,s0) 
A„  S0 


"2  lnCX(j)]  *  X2n 


Decision  Rule:  Decide  N  +  S  iff 

T>4,l-<»/2  0r  T<4,c./2 
Test  the  following  PN  situation  using  generated  Pareto  data: 

PN:  (A,s)  =  (1,2)  vs.  N  +  S:  (A,s)  t  (1,2) 

Data  Set:  (see  Table  2) 

1.  Xj,  ...»  X50  i.i.d.  Pa(l;2) 


Decide  N  +  S  iff  T  >  140.2  or  T  <  82.4 

Since  the  calculated  value  of  the  test  statistic  is  T  *  91.02,  one 
decides  PN. 

Example  6.6.  A  unknown  s  unknown,  /.( A,s)  e  ft(PRP) 

PN:  s  <_  Sp  vs.  N  ♦  s :  s  >  sg 

n  2 

Detector  Statistic:  T  ■  2sQ[  \  In  X.  -  n  In  X(l)]  'v  X2n_2  Q 

j-1  J  ’ 

Decision  Rule:  Decide  N  ♦  S  iff  T  <  xiL  -  ~ 

—  - - ■  2TI-2,(X 
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I 

i 


4 


4 


I 


Test  the  following  PN  situation  using  generated  Pareto  data: 

PN:  s  <  3  vs.  N  +  S:  s  >  3  (a  =  .05) 

Data  Sets :  (see  Tables  1  and  2) 

1.  Xj,  ....  Xg0  1.1. d.  Pa(l;2) 

2.  Yj,  ....  Y50  i.i.d.  Pa(5;5) 

2 

The  critical  value  is  X^g  =  76.5 
The  calculated  statistic  values  are 

From  data  set  1:  T  *  129.32,  so  one  decides  PN. 

From  data  set  2:  T  *  72.64,  so  one  decides  N  +  S. 


Example  6.7.  A  unknown,  s  unknown  L(A,s)  £  ft(PRP) 

PN:  A  <  Aq  .  vs .  N  ♦  S:  A  >  AQ 

n(n-l) [In  X(l)  -  In  A  ] 

Detector  Statistic:  T  «  -  ^  F,_ 

S 

2,  In  X.  -  n  In  X(l) 

1  5 

Decision  Rule:  Decide  N  +  S  iff 

i.  X(l)  >  Aq  and 

ii.  T  >  f ,,  - 

(2,2n-2,l-a) 

Test  the  following  PN  situation  using  generated  Pareto  data: 


,2n-2) 


PN:  A  <  3 


vs .  N  S :  A  >  3 


Data  Sets :  (see  Tables  1  and  2) 

1.  Xj,  ...  Xg0  l.i.d.  Pa(l;2) 

2.  Yj,  ....  Yso  i.i.d.  Pa(5;5) 


The  critical  value  for  T  (a  =  .01)  is 


f(2,98, .99)  =  4-87 


Since  X(l)  *  1.024,  one  decides  PN  for  data  set  1. 

Since  Y(l)  *  5.006,  and  Ty  *  103.7,  one  decides  N  +  S  for  Data 
Set  2. 


>.8.  A,  s  unknown 


PN:  L(A,s)  e  Jl(PRP)  vs.  N  ♦  S:  L(A,s)  $  Sl(PRP) 


Detector  Statistics: 


i  i  n_ *■ 

i.  D*  »  sup  l  e  (z  -  pJ — )  -  z|  'v  K-S(n-2) 


2  l  n"2  E-i 

ii.  D^  »  sup |—— «•  l  t  (z  -  ^ — )  -  z|  'v  K-S(n  -  2) 
n  z  n'£  j-l  cn-l 


where 


•  1  *"  »  *  <»  • 


1  <  r  <  n  -  1 


K  •  Z  (n  ♦  1  -  j) 


1  <  r  <  n  -  1 
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Decision  Rules:  Decide  N  ♦  S  iff 


i.  D  >  d,  «  .. 

n  (n-2,ot) 


ii.  D  >  d,  -  . 

n  (n-2,a) 


where  d(n  -  2, a)  is  appropriate  value  from  Kolmogorov- Smirnov  table. 

Data  Sets:  (see  Tables  1,  2  and  3) 

1.  Xj,  ....  X5Q  i.i.d.  Pa(l;2) 

2.  Yj,  ....  YS(J  i.i.d.  Pa(S;S) 

3.  Wj,  ....  W24  waiting  times  for  "peaks"  (i.e.  R's)  of  heart-beat 

cycles . 

The  critical  value  for  data  sets  1  and  2  is: 


d(48, .01)  *  0,22 


The  critical  value  for  data  set  3  is: 


d(22,.01)  '  °-314 


The  calculated  statistic  values  are: 


“so  <*>  ■  °-058 

Decide 

PN. 

°S0<«  -  0.058 

Decide 

PN. 

Dcn(Y)  «  0.042 
*>u  ^ 

Decide 

PN. 

080(Y)  .  0.043 

Decide 

PN. 

o‘4(ir)  -  0.614 

Decide 

N+S 

D2  (W)  ■  0.605 

*4  *\j 

Decide 

N+S 
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Example  7.1.  Aj,  known;  s^,  unknown 


PN:  s  *  s 2  vs .  N  +  S;  s ^  t  s ^ 


m 

n  l  ln(X./A  ) 
i  =1  ^  1 

Detector  Statistic:  T  *  - ^ - -  'v  F(2m,2n) 

m  l  ln(Y  /A  ) 

j-1  3  * 


Decision  Rule:  Decide  N  +  S  iff 


T  >  F(2m,2n,l-a/2)  or  T  <  F(2m,2n,a/2) 


Data  Sets:  (See  Tables  .  1  and  2) 


1. 

V 

....  x50 

i.i.d. 

Pa(l;2) 

2. 

Yr 

Y 

x50 

i.i.d. 

Pa(5;2) 

3. 

zr 

Z50 

i.i.d. 

Pa(5;5) 

The  critical  values  are: 

f  (100, 100 ,.995)  1,68 
f(100,100, .005)  =  0,595 

The  calculated  statistic  values  for  each  pair  of  data  sets  are  given  below 


From  data  sets 

1 

and 

2: 

T  =  0.788, 

so 

one 

decides 

PN. 

From  data  sets 

1 

and 

3: 

T  =  1.868, 

so 

one 

decides 

N  ♦  S. 

From  data  sets 

2 

and 

3: 

T  «  0.422, 

so 

one 

decides 

N  ♦  S. 
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Example  7.2.  A^,  unknown;  s^,  known 


PN:  Al  =  A2 


vs.  N  t  S:  Aj  +  A2 


Detector  Statistic:  T  =  In  (tt£;4-)  ^  F(x) 

A 


where 


F(x) 


ms  -ns.x 

— - e  x  >  0 


(msj+ns^ 


ns-  ms  x 

e  x  <  0 


(ms1+ns2) 


Decision  Rule:  Decide  N  +  S  iff 


>  Cj  or  T  <  C2  where  FfCj)  *  1  -  a/2 


F(c2)  =•  a/2 


Data  Sets :  (see  Tables  1  and  2) 


1. 

V 

....  X50 

i.i.d.  P*(l;2) 

2. 

V 

Y 

**  »  50 

i.i.d.  Pa(5;2) 

3. 

V 

Z50 

i.i.d.  Pa(5;S) 

The  first  two  data  sets  constitute  a  special  case  (Sj  *  $2  *  s0  known) 
and  will  be  treated  in  Example  7.2a.  For  the  remaining  two  pairs  of 
data  sets,  one  has  the  following  critical  values:  (m  ■  n  *  50, 

Sj  *  2,  s2  =  5,  a  ■  .01) 

Cj  *  0.016  and  c2  »  -0.050 


The  appropriate  decision  rule  is: 


Decide  N  +  S  iff  T  >  0.016  or  T  <  -0.050. 


The  calculated  statistic  values  are  given  below: 

From  data  sets  1  and  3:  T  =  =  1-587,  so  one  decides  N  +  S 

From  data  sets  2  and  3:  T  *  In (^^~)  =  -0.008,  so  one  decides  PN. 


Example  7.2a.  A^,  unknown;  si  =  s2  =  s0  known 


PN:  Aj  =  A2 


vs.  N  ♦  S:  Aj  f  A2 


Detector  Statistic: 


t.IQI 

T  X(l) 


Decision  Rule:  Decide  N  *  S  iff  T  <  bj  or  T  >  b2 

u  „  l/msn  M  „  -l/ns_ 

where  bj  ■  £  (§)]  ,  b,  -  §  <§)] 


Data  Sets:  (See  Table  1) 

1.  Xj,  ....  X5Q  i.i.d.  Pa(l;2) 

2.  Yj,  ....  Ysq  i.i.d.  Pa(5;2) 

The  critical  values,  for  o=.01  are: 

bj  ■  0.95  and  b2  *  1.05 

The  calculated  statistic  value  is 

T  =  -  4.929  Decide  N  +  S. 


J 


! 

I 

I 


J 


1 
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Example  7.3.  =  A  unknown,  s^,  S£  unknown 


PN: 


S1  —  S2 


Detector  Statistic: 


vs .  N  +  S:  3^82 


where 


m 

n  =  1  In  X.  -  a  In  X(l) 
j=l  3 


and 


m  n 

p  ■  l  In  X.  +  l  In  Y.  -  NW 

j*l  3  j»l  3 


where  W  *  min  (in  X(l) ,  In  Y(l)> 

Decision  Rule:  Decide  N  +  S  iff  T  <  c  where  c  is 


Ic(m  -  1,  n)  *  1  -  a 

Data  Sets :  (See  Tables  1  and  2) 

1.  Xj,  Xg0  i.i.d.  Pa(5;2) 

2.  Yj,  ...,  Yg0  i.i.d.  Pa(5;5) 

Computation  of  critical  value: 

I  (49,50)  »  0.99  (take  a  *  .01) 

c 

c  *  0.61  (from  tables  of  Incomplete  Beta  function) 


The  decision  rule  is: 

Decide  N+S  iff  T  <  0.61 


.694 


I 


determined  by 
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one  decides  PN. 


Example  7.4.  A^  A 2  unknown;  Sj  ■  s2  =  s  unknown 


PN:  Ax  £  A2 


vs •  N  +  S:  A^  >  A^ 


Detector  Statistic:  T  =  In  X(l) 

Decision  Rule:  Decide  N  +  S  iff  T  >  c  where  c  = 


nc 

n 


W  =  min  {In  X(l) ,  In  Y(l)> 


H  *  [  In  X.  ♦  l  In  Y. 
j*l  3  j»l  3 


-  NW 


and  c'  is  determined  by 


(J)  (1  -  c')N"2  -  a 


Data  Sets :  (See  Tables  1  and  2) 

1.  Xj,  ....  XSQ  l.i. d.  Pa(l;2) 

2.  Yj,  ....  Y$0  i.i.d.  Pa(5;5) 


Computation  of  critical  value  (a  «  .01): 


«TO)C1  -  c')98 


.01 


.04 


since  h  »  129.73  and  W  ■  In  X(l)  *  .024,  one  has 


c  .  +  ,024  «  .128 
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since  T  *  .024,  one  decides  PN. 

Example  7.5.  Aj,  A2  unknown;  Sj, 

PN:  s2  “  s2  vs.  N  +  S: 

(n  -  2)n 

Detector  Statistic:  T  ■  - jyjy 


s 


1 

2 


unknown 


%  F  (2(»-l),  2(n- 


m-1 


"here  ni  ‘  jl,  ln 


and 


„  .  V  in  Xii  *  » 

n2  fa  111  ni) 


Decision  Rule:  Decide  N  ♦  S  iff 

T  >  f(2(m-l) ,2(n-l) ,  (l-a/2}) 
or 

T  <  f(2(*-l),2(n-l),tf2) 

Data  Sets:  (see  Tables  1  and  2) 

1.  Xj ,  ....  Xjq  i.i.d.  Pa(l;2) 

2.  Yj ,  ....  Y50  i.i.d.  Pa(5;2) 

3.  Zj,  ...,  Z5Q  i.i.d.  Pa(5;5) 

Critical  values:. 

^(98,98,0.99)  -  1,69 
f (98,98,0.01)  “  0,592 


\ 


B 


- 

* 


i 


t 


> 

« 

r 


> 

. 


i 


r 


» 
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The  decision  rule  is: 


Decide  N  ♦  S 

iff 

T  >  1.69 

or  T 

<  .592 

For  data  sets 

1  and 

2: 

T  • 

0.760, 

so 

one 

decides 

PN. 

For  data  sets 

1  and 

3: 

T  - 

1.779, 

so 

one 

decides 

N  ♦  S. 

For  data  sets 

2  and 

3: 

T  » 

0.427, 

so 

one 

decides 

N  ♦  S. 

Example  7.6.  No  techniques  for  Case  7.6  have  been  developed  by 
the  authors. 

Example  7.7.  Aj,  Aj  unknown;  sj*  »2  unknown 

PN:  (A1,  Sj)  -(Aj,  s2)  vs.  N  ♦  S:  (Aj,  Sj)  ^  (A^  s2) 

Detector  Statistics 
(n  -  2)n. 

(i)  ?1  -  (>  _  ^  *  p(2(«-l),2(n-l)) 

<■  -  -  ti«  t2  .  anfcam«H' 

Decision  Rule:  Decide  N  ♦  S  iff 

(i)  T1  >  f(2(«-l),2(n-l),l-a/4)  or 
T1  <  f(2(»-l).2(n-l),a/4) 

(ii)  If  fa/4  <  T,  <  tlma/A 

decide  N  ♦  S  iff 

T2  >  f(2,2N-4,l-o/2) 


1 


Data  Sets :  (see  Tables  1  and  2) 


1. 

V 

••••  X5o 

i.i.d. 

Pa(l;2) 

2. 

V 

Y50 

i.i.d. 

Pa(5;2) 

3. 

V 

Z50 

i.i.d. 

Pa(5;5) 

Critical  values  (a  =  .01): 

Cl)  f (98, 98,. 9975)  “  1,77 
f (98, 98,. 0025)  *  *S6S 

(ll)  f (2,196, .995)  "  5,35 
Decide  N  +  S  iff 

(i)  Tj  >  1.77  or  T  <  .565 

(ii)  If  .565  <  T1  <  1.77,  decide  N  ♦  S  iff  T2  >  5.35 


1. 

For 

ex, 

'X/ 

Y). 

*\j 

T1 

■*  0.760, 

V 

■  39.102  so  one  decides 

N  + 

2. 

For 

(X, 

'V/ 

Z), 

T1 

»  1.779, 

T2 

57.682  so  one  decides 

N  ♦ 

3. 

For 

Of. 

% 

Z), 

T1 

*  0.427, 

T2  * 

0.251  so  one  decides 

N  +  S 

sDuma  am  iznaa 
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rnmmHtmmt) 


m 

i 


to.  «(V  *0*00  (*■ 


i  tf  «r  ««• 


Pareto  distribution,  renewal  process,  distribution-free;  sufficient 
statistic;  statistical  noise;  order  statistics;  Kolmogorov -Smirnov 
statistic;  conditional  distribution;  unbiased,  Lilliefor's  statistic, 
Srinivasan  statistic. 


ami  a  act  (Crnmmmm  ****  mtm  it  wihim  «*o  tmmtiir  tf  nma 

7  Employing  minimal  sufficient  statistics  maximal  statistical  noise, 

'  several  Kolmogorov-type  statistics;  and  conditional  distributions, 
optimal  detection  procedures  are  constructed  for  various  one-and-two- 
sample  problems  Involving  Pareto  Renewal  processes.  Cases  with  and 
without  nuisance  parameters  are  treated.  Optimal  parametric  and 
distribution-free  procedures  are  developed,  r- — - 
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